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ABSTRACT 

We present an analysis of the science image quality obtained on the twin 6.5 m Magellan telescopes 
over a 1.5 year period, using images of ^ 10^ stars. We find that the telescopes generally obtain 
significantly better image quality than the DIMM-measured seeing. This is qualitatively consistent 
with expectations for large telescopes, where the wavefront outer scale of the turbulence spectrum 
plays a significant role. However, the dominant effect is found to be wind speed with Magellan 
outperforming the DIMMs most markedly when the wind is strongest. Excluding data taken during 
strong wind conditions (> 10 m s~^), we find that the Magellan telescopes still significantly outperform 
the DIMM seeing, and we estimate the site to have £o '^ 25 m on average. We also report on the first 
detection of a negative bias in DIMM data. This is found to occur, as predicted, when the DIMM is 
affected by certain optical aberrations and the turbulence profile is dominated by the upper layers of 
the atmosphere. 

Subject headings: Astronomical Phenomena and Seeing - Astronomical Instrumentation - Astronom- 
ical Techniques - Data Analysis and Techniques 
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1. INTRODUCTION 

Much time and money is invested in finding optimal 
sites for the next generation of big optical-IR telescopes 
and it is of interest to learn the detailed atmospheric 
properties of major observatory sites since the seeing is 
a major limitation to the resolution and signal to noise 
ratio that can be achieved. It is well known that large 
aperture telescopes frequently perform better than we 
would naively expect from the simple seeing measure- 
ments at ma jor obser vatories (e.g. Mart i n et al. 1998.; 
Wilson et ai]|l 999: To kovinin et all 120071: ISalmon et atl 
20091: lEls etal., ,2009) in terms of image full width at 
half maximum (FWHM). Most astronomers consider 
the "seeing" to be synonymous with the FWHM of 
a stellar image on their detector, but this assessment 
is formally incorrect. Seein g is a sim ple measure- 
ment of Fried's parameter, tq (jFriedl [19651 ) . determined 
usi ng differential image mot ion monitors (DIMMs - 
see lSarazin and Roddieijri990l and references therein) but 
there are other important parameters governing imaging 
quality, such as the characteristic time constant, the iso- 
planatic angle and the wavefront outer scale, Cq. These 
can play important roles in different regimes. 

1.1. Background: Turbulence, seeing and outer scale 

Using a single parameter to determine image FWHM 
implicitly assumes that energy is injected into the atmo- 
sphere at infi nite scale following the Kolmogoro v turbu- 
lence model (jTatarskiil 119691 : lKolmogorovlll941[ ). an ap- 
proximation that is valid only for small telescopes. In 
the Kolmogorov model, energy is gradually transferred 
to smaller and smaller scales until it is dissipated through 
the viscosity of the air on scales of ~ mm - cm (the inner 
scale, Iq). The result is local variations in atmospheric 
temperature and density, and thus in refractive index, 
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n. In reality the atmosphere continually has energy in- 
jected into it through heating, convection currents and 
wind shear on scales Lq ~ 10 m, which we can introduce 
as a cutoff in the turbulence spectrum or "outer scale 
of turbulence" , Cq . This is known as the von Karma n 
tmbulence model (see lBorgninolll990t IZiad et al.ll2000( ). 
In practice atmospheric turbulence results in small 
amounts of image motion. The standard device for mea- 
suring this motion is the DIMM. A DIMM measures the 
variance of the differential image motion, a'^, of two im- 
ages of the same star, observed through separate aper- 
tures in the entrance pupil of a small telescope. This 
gives a measurement of tq for a sub-aperture of diam- 
eter D and wavelength A: cr-^ ex Xr^ D~^/^. Thus 
the FWHM eq of the long-exposure seeing-limited point 
spread function in large telescopes is computed within 
the Kolmogorov turbulence model as eo = 0.98A/ro (x 

(i?/A)V5(,2)3/5. 

From a turbulence spectrum it is possible to derive 
a statistical structure function F^(r) which describes 
the phase decorrelation between any pair of points in 
the atmosphere, separated by r. In the Kolmogorov 
model, the phase structure f unction (the in ean square 
differential phase error - see iRoddien Il981l for a use- 
ful derivation) grows larger with distance of separation, 
Dcf,{r) ex (|r|/ro)^/^ - see dotted line in Fig. [H But 
in the von Karman mod el the structure function be- 
comes more complicated (jConsortini and Ronchil |1972|) 
and grows rapidly with r up to r = Cq where it saturates 
- see solid line in Fig.[T] Beyond Co, D^{r) no longer in- 
creases. The outer scale effect is detectable in telescopes 
that are an appreciable fraction of the size of Cq , where 
it begins to deviate from the Kolmogorov case. 

From the structure function we can derive the opti- 
cal transfer function (OTF) of a large telescope, To(/) = 
exp[— O.5-D0(A/)], where /is the spatial frequency in in- 
verse radians. The corresponding Point Spread Function 
(PSF) is found by taki ng the Fourier Transform of the 
OTF. iTokovininI (|2002t ) obtained the PSF for the von 
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Fig. 1. — The von Karinan phase structure function, D^{r) (soUd 
hne Co = 25 m, dashed line Co = 1000 m) compared to the Kol- 
mogorov phase structure function (dotted hne), aU for ro = 0.2 m. 
Unhke the Kolmogorov structure function which continues to grow 
with increasing r, the von Karman structure function saturates at 
r =Co- The outer scale effect becomes noticeable at scales signifi- 
cantly smaller than Co- 



Karman case by renormalizing the OTF. The resuh is 
an approximation of the seeing-limited PSF of a large 
telescope, with FWHM c^k approximated by: 
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which we use later in this paper to estimate an average 
value of £o at the Las Campanas site. The effect of 
introducing an outer scale to the phase structure function 
is to decrease the FWHM of the point spread function 
(PSF) of a "large" telescope. 

2. IMAGE QUALITY MEASUREMENTS AT THE 
MAGELLAN TELESCOPES 

The twin 6.5 m Magellan Telescopes (individually 
named Clay and Baade) were built and are operated 
by a consortium consisting of the Carnegie Institution 
of Washington, Harvard University, MIT, the University 
of Michigan, and the University of Arizona. They are 
situated on Cerro Manqui within the Las Campanas Ob- 
servatory site at (29°00'54"S 70°41'32"W) on the south- 
ern edge of the Atacama desert in Chile. Cerro Man- 
qui sits at 2,450 m, slightly lower than Las Campanas 
peak itself at 2,551 m. The telescope s are described 
in detail by 'Shect man and Johns! (|2003[ ) , and the Mag- 
ellan optical design including its active opt ics Shack- 
Hartmann (S-H) system is further discussed bv Sh ectmanI 
(|1994D . Optic al seeing down to a FWHM ~ 0'.'25 have 
been recorded (jShectman and Johnsll2003l ). 

In 2007 we began monitoring the image quality on all 
Magellan telescope science images with > 30 s exposure 
time. The aim was to see how well they do compared 
to the nominal seeing conditions and to help diagnose 
misalignments in the optical path. Furthermore, Las 
Campanas Observatory (LCO) is the chosen location for 
the Giant Magella n Telescope (GMT) for which site test- 
ing b egan in 2005 (|Consortiumll2006l : lThomas-Osip et al.l 
I2008f) . The site testing for GMT provides us with a 
wealth of tools for determining the precise seeing, atmo- 
spheric turbulence and meteorological conditions - see 
section 12.21 The GMT site testing program made use 
of four DIMMs around the Las Campanas site including 
one for reference at the Magellan Telescope site on Cerro 
Manqui. 



In this paper we take the image-wide average full 
width at half maximum (FWHM) of Magellan science- 
camera stellar images using a program called "DIPSF" 
(as described in section [2. II) and compare them with the 
contemporaneous DIMM (see section [2?2|) and Magellan 
guide camera FWHM measurements of the seeing. We 
then consider simultaneously taken atmospheric turbu- 
lence and meteorological data, using the GMT Multi- 
Aperture Scintillation Sensor (MASS) and weather sta- 
tion. 

The Magellan guide cameras use an RG610 filter, with 
a mean wavelength of 765.2 nm. DIMM measurements 
are calculated at 500 nm. The guide camera and DIPSF 
FWHM measurements in this paper are corrected to air- 
mass 1 and to a wavelength of 0.5 /xm following the well 
know To oc A^/^ (eo oc A~-^/^) and cq oc airmass'^^'^ re- 
lations. The guide camera measurements do not neces- 
sarily refiect the actual seeing at the telescopes since the 
guider optics are not capable of producing images better 
than 0'.'4 arcsec (FWHM). The guide stars are several 
arc minutes off-axis and are uncorrected for aberrations 
in the telescope optics. 

2.1. Image Analysis Software 

We use a source detection and fitting program known 
as D IPSF, which is based on DoPhot (Schechter et al.l 
Il993f ). DIPSF reads in an image and bitmaps it to 
form a low dynamic range image in memory, which is 
then analyzed for sources. Only the most basic data 
reduction steps are performed on the image: Bias sub- 
traction is performed using the overscan regions, and 
static bad pixel masks (computed during engineering 
runs each month) are applied, including any vignetted 
regions of the detector. We stress that the raw data 
remains untouched - the data is read in and all process- 
ing is performed in memory. Sources are are fit with a 
2-dimensional gaussian and are classified as stellar, cos- 
mic ray, double or extended according to the resulting fit 
parameters. The contemporaneous guide camera image 
FWHM is used as a guess for the first fit, and for ref- 
erence in decision making about the nature of a source. 
Each star detected on a given frame is recorded in an 
output file, and several high level statistical measures of 
the average point spread function across the image are 
stored to a database. 

Magellan data are proprietary and are stored for just 
one month in case of corruption of a dataset that an 
astronomer takes home. Therefore each image is ana- 
lyzed the day after it was taken. At present DIPSF is 
run by the instrument specialists on each data directory 
from a night's observing. Only direct images that have 
exposure times > 30 s are processed in order to ana- 
lyze seeing-averaged images. We run DIPSF on data 
from LDSS3 (on Clay ) MagIC and IMACS (on Baade) 
- see 10 sip et al.l ()2008l 120041) for more information on the 
instruments themselves. There is also a mode for reading 
MegaCam data. 

The code is written in Fortran and is called from a 
Python script. Output includes list of all sources identi- 
fied (-STAR files), along with their properties: position, 
size, ellipticity, position angle, flux. Two high level out- 
put files contain the average PSF characteristics over 
each chip (for multi-chip instruments) and over the en- 
tire frame (_CHIP and _frame files respectively). Af- 
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Fig. 2. — Example DIPSF output showing all stars detected on a 
crowded IMACS f/2 image from the night of 2009 April 23rd. The 
ellipses denote each stellar PSF and are shown with true ellipticity 
and position angle. Ellipse size is scaled according to the FWHM 
of the stellar image, and each ellipse is color-coded according to its 
relative size. On this occasion the mean image FWHM is 0'.'504 ± 
0''014. Axis labels show pixel number measured across the IMACS 
mosaic. 



ter running, the python script posts the contents of the 
frame-averaged results to an SQL database. This can 
then be queried along with the meteorological and guide 
camera data from any given night or nights. The indi- 
vidual -STAR output files are stored on a server. The 
contents of an example IMACS f/2 _STAR file are plotted 
in Fig. [2] to highlight the stability of the PSF across the 
entire half-degree field of view. Stellar images are shown 
to scale with their true ellipticities and position angles 
and are color-coded according to FWHM. In this image 
the mean image FWHM is 0'.'504 ± 0'.'014. The software 
can also be run in a "diagnostic" mode which runs on 
all images, the results of which are simply stored locally, 
not posted to the database nor server. This is useful for 
analysis of engineering data, such as pinhole masks and 
out of focus "donut" images, and may be used to detect 
and diagnose defects in the telescope optical path, and 
potentially to supply Shack-Hartmann like information 
to the telescope mirror actuators. 

The program was run from May 2007 on occasional ba- 
sis, and on all data from the start of 2008 to the present. 
We thus have a database containing the frame-averaged 
PSF FWHM for every > 30 s LDSS3, IMACS and MagIC 
image taken on the Magellan telescopes since January 
2008, and a few images from earlier. The software con- 
tinues to be run on a regular basis, but the DIMM's are 
now infrequently operated since the GMT site testing 
program is complete. 

2.2. GMT site testing: DIMM seeing, meteorology and 
turbulence measurements 

An extensive site testing program has identi- 
fied Cerro Las Campanas as the future GMT 
site (jThomas-Osip et al.l I2010D . Four sites within the 
LCO property were studied. Of these, the site of the 
Magellan telescopes, Cerro Manqui, was used for ref- 



erence purposes only. The seeing at all four sites was 
monitored through the use of identical DIMM systems 
in 7 m towers. The DIMM method is based on the re- 
lationship between the FWHM of a long exposure image 
in a large telescope to the variances in the differences in 
the motion of two images of the same star through the 
use of Kolmogorov theory. The two images are created by 
placing a mask with two sub-apertures containing prisms 
at the entrance of the optical tube. The GMT instru- 
ments are based on the CTIO RoboDIMM^ but have 
several improvements. Image quality has been improved 
by using two thinner prisms as opposed to one thicker 
prism and an open aperture. Following a technique de- 
veloped by the Thirty Metre Telescope (TMT) project, 
the Carnegie DIMM (CDIMM) software uses a drift scan 
readout mode which allows for many more image motion 
measurements to be made per minute and thus improved 
statistics^. The Manqui DIMM is located approximately 
20 m northwest of the Clay dome, almost exactly abeam 
of the Magellan telescopes with respect to the prevailing 
northeasterly wind. The Manqui peak itself is approxi- 
mately 3.8 km at bearing 345° from Campanas peak. 

The DIMM measurements are filtered to remove poor 
quality data using three criteria: focus, read noise, and 
the number of measurements taken during one minute. 
We measure focus using the mean separation of the two 
images during one minute of measurements. Nominal fo- 
cus is determined on a night with good seeing by varying 
the focus and maximizing the Strehl ratio. The CDIMM 
software is designed to keep the mean separation within 
a range of 1 pixel on either side of the nominal focus. 
There are however some measurements throughout the 
night (but especially at the beginning of the night) that 
need to be removed because the separation is outside 
the nominal range for a short time while the instrument 
works to correct it. The filters for read noise and the 
number of measurements allow for the removal of sus- 
pect quality data due to clouds or tracking errors. 

In addition to seeing, many other parameters were 
also studied as a part of the GMT site testing program. 
Some of these will also be used here in order to under- 
stand differences seen between the Magellan telescope 
(DIPSF) FWHM and the DIMM seeing FWHM. Me- 
teorological characteristics (pressure, temperature, hu- 
midity, and wind) were collected using weather stations 
manufactured by Davis Instruments Corp and mounted 
on 10 m towers. Turbulence in the at mosphere above 
500 m is being monitored by a MASS (jKornilov et al.l 
l2003f ). The spatial scale of the scintillation variation 
depends on the distance to the layer in which the tur- 
bulence giving rise to the wave front phase disturbance 
exists. Thus, the turbulence profile, denoted C^{z) dz, 
at a small number of discrete layers (z = 0.5, 1, 2, 4, 8, 
and 16 km) can be restored by fitting a model to the dif- 
ferences between the scintillation indices within four con- 
centric apertures. In addition to the turbulence profile, 
the MASS also measures free atmospheric seeing (essen- 
tially the integral of the turbulence profile). The differ- 
ence between the total DIMM and the free atmosphere 
MASS turbulent integrals is a measure of the portion of 
the total seeing contributed by a ground layer below 500 
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m (jTokovinin and Korni"Iovll2007[ ). 

The MASS data have been reprocessed usin g a new 
version of Turbina (jKornilov and Shatskvl 120051 ) to cor- 
rect for the effect of strong s cintillation under poor seeing 
conditions (Tokov ininI l2006l ) . Strong scintillation affects 
the turbulence profile and free atmosphere seeing by over- 
estimating turbulence and spreading turbulence to lower 
levels. The sensitivity of the possible errors due to in- 
correctly determined parameters (such as the instrument 
magnification or the non-Poisson and non-linearity char- 
acteristics of the d etect ors) lead to a maxim um bias on 
the order ±0'/ 05 (iTho mas-Osi p et all 120081 ) similar to 
that found bv lEls et~ai.i (|2008i) . The seeing, e, contri- 
bution by the ground layer (defined here as lower than 
500 m) can be inferred from the MASS and DIMM mea- 
surements in the following manner: 
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The absolute value is used only in computing statisti- 
cal distributions. This need arises due to that fact that 
some small percentage (2-20%) of matched MASS and 
DIMM measurements result in what could be termed a 
negative ground layer seeing values in the sense that the 
MASS measurement is larger than the DIMM measure- 
ment. This is, of course, not physically possible but is a 
result of non-zero errors and the subtraction of quanti- 
ties of similar magnitudes. The range in the percentage 
of negative ground layer seeing arises when considering 
selections of data filtered for scintillation index in the 
smallest MASS aperture. The correction made in the 
MASS reprocessing for strong scintillation is not perfect 
and this effect can be seen by examining the percentage 
of negative ground layer found in samples with varying 
scintillation indices. The sample presented in this anal- 
ysis is not filtered for scintillation index since this would 
introduce a bias towards better seeing. But in this case, 
one must realize that the ground layer as determined 
from data that has not been filtered for scintillation in- 
dex will be less accurate. 

2.3. DIMM bias 

We identified a problem in the DIMM data, start- 
ing in January 2009. The Meade telescope primary 
mirror developed a non-repeatable movement during fo- 
cus changes which caused the collimation to deteriorate. 
This was discovered in the consequential deterioration of 
the DIMM Strehl ratios. We therefore restrict this pa- 
per to observations taken between July 2007 and the end 
of December 2008. The bias and its effects on the data 
taken in 2009 are discussed further in section [ 



3. THE DIPSF DATASET AND MATCHING TO THE DIMMS 

We have cut out the 2009 data due to the DIMM prob- 
lem mentioned in section [2731 DIPSF data was obtained 
on 230 individual nights over the period July 2007 - De- 
cember 2008. We have identified 70 nights during this 
period on which there were problems with the DIPSF 
data, due to source detection or other software errors, 
or out-of-focus observations. Fits to '^ 7 x 10^ stellar 
images were performed during this time, yielding 5617 
frame-averaged FWHM datapoints with an overall av- 
erage of 128 stars per data point. The distribution by 
instrument is summarized in table [T] For all data we 



TABLE 1 
Summary of results: Number of fields, np, analyzed for 

EACH INSTRUMENT, AND AVERAGE NUMBER OF STARS PER FIELD, 

rf*. The final column lists the number of fields with 

CONTEMPORANEOUS DIMM DATA, ricD ■ 



Instrument Field area (sq. arcniin) 



np 



n^o 



IMACS 
Magic 
LDSS3 



635 


1287 


459 


199 


6 


2934 


17 


310 


54 


1396 


42 


367 



have a corresponding guide camera image FWHM value, 
as this is required input for the DIPSF code (see sec- 
tion [2T] above). In order to be useful for further analysis 
we also require a contemporaneous DIMM measurement. 
We query the DIMM database once for each night on 
which we have DIPSF data. Once we have downloaded 
the entire night's DIMM data it is straightforward to 
match them to our DIPSF exposures (and much faster 
than performing each individual match in SQL directly 
on the database). We average all DIMM readings taken 
during each DIPSF exposure to give the corresponding 
DIMM seeing for that exposure. An identical approach 
is taken with the guide camera readings. Unfortunately, 
the Cerro Manqui DIMM was not operated every night 
during the time period for which the DIPSF data was col- 
lected, as some of this period occurred after the GMT site 
testing program finished its main data collection. There- 
fore we do not have corresponding DIMM seeing values 
for the majority of DIPSF data points. After match- 
ing we have 876 concurrent data points from 62 nights 
(^ lO^stars). The breakdown by instrument is shown in 
table [1] Fig. [3l shows examples of the variation of image 
quality throughout various nights. The DIPSF science 
image quality can be seen to track the DIMM and guide 
camera seeing values very well. 

4. RESULTS 

In Fig. m we plot the DIPSF seeing values against 
simultaneous readings for the DIMM seeing (left) and 
guide camera seeing (right). It is clear that the DIPSF 
and DIMM values are strongly correlated (Spearman's 
rank correlation coefficient p = 0.71; p = 0.00). There 
is however a large dispersion, and a deviation from 1:1 
(see left panel of the figure). The DIPSF and guide cam- 
era image qualities are also strongly correlated, and lie 
closer to 1:1, although with a larger dispersion (see right 
panel of the figure). The DIPSF and DIMM FWHM dis- 
tributions are significantly different according to a two- 
sided Kolmogorov-Smirnov (K-S) test, which returns a 
cumulative difference D = 0.21 between the two distri- 
butions, with a probability p ^ 10^^^ of being drawn 
from the same underlying distribution. The difference is 
even more significant for DIPSF and guide camera image 
quality distributions, with a K-S test returning a cumu- 
lative difference difference D = 0.33 with a probability 
p~ 10-42. 

The median DIMM seeing for our matched dataset is 
0'.'625. The median DIPSF FWHM is 0'.'575 and for the 
guide cameras it is 0'.'688. Mean values and standard 
deviations are FWHMdimm = 0'.'658 ± 0'.'205 (skewness 
= 1.04), FWHMdipsf = 0'.'585±0'.'150 (skewness =1.24) 
and FWHMguide = 0'.'697 ± 0'.'165 (skewness =0.923). 

The Magellan image quality is better than the DIMM 
seeing 69% of the time, and better than the guide cam- 
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DIPSF, DIMM and guidecam seeing: 2454417 
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Fig. 3. — Examples of contemporaneous DIPSF (black squares), DIMM (dark grey crosses) and guide camera (light grey crosses) seeing 
measurements for (ut-dates) 2007-11-12, 2007-11-14, 2007-11-15 and 2008-01-16. 



era image quality 90% of the time. In poor seeing 
conditions (DIMM values > 1") Magellan does better 
than the DIMM 98% of the time. The result is qual- 
itatively as we would expect for a large telescope: In 
general the Magellan science cameras return an image 
quality that is better than the seeing as measured by 
the DIMMs. At good seeing we are limited by the 
telescope optics, and the two readings converge, while 
at poor seeing, the Magellan telescopes appear to do 
progressively better than would be expected from the 
DIMM seeing. The best linear fit to the entire dataset 
gives FWHMdipsf = 0.53 FWHMdimm + 0.24. For the 
individual instruments this varies little (see Fig. HI left 
panel), with LDSS3 following a slightly steeper relation 
than Magic and IMACS. When comparing the DIPSF 
data to the guide camera FWHM there is a more marked 
difference between the instruments. 

It is interesting to note that the IMACS DIPSF 
FWHM values do not get as low as those for the other in- 



struments, in spite of having the largest number of mea- 
surements. The minimum IMACS DIPSF FWHM (ex- 
cluding a single outlier at 0'.'24) is at 0'.'35. LDSS3 has a 
minimum FWHM at 0'.'22 and 17 points that are lower 
than the IMACS minimum (excluding a single outlier 
at O'.'IT). Magic has a minimum FWHM of 0'.'27 and 
13 points that are below the IMACS minimum. This 
may be traceable to the averaging over a larger field 
of view, or inherent to the optics on the IMACS port 
of the telescope. The (S-H) wavefront sensors fitted to 
the instrument ports provide a measure of the optical 
aberations present during operation of the Magellan tele- 
scopes. These predict a modal Gaussian PSF FWHM 
of O'.'IT. This is sufhcient to degrade a seeing of 0'.'25 
to O'.'SO. However, the S-H system on IMACS predict 
a modal Gaussian PSF FWHM of 0'.'30 FHWM. This is 
sufficient to degrade 0'.'25 seeing to 0'.'39. In Fig.[5]we plot 
the cumulative histograms the DIPSF, DIMM and guide 
camera FWHMs. On the right hand figure, we show the 
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Fig. 4. — FWHM as measured by DIPSF on science images against simultaneous DIMM seeing FWHM values (left) and guide camera 
seeing FWHM values (right). Red crosses denote IMACS data, green squares LDSS3 and blue triangles MaglC. The solid line on each 
panel shows 1:1 for reference. We also show best linear fits for the whole sample (black dashed line) and for each instrument individually 
(dotted lines in same colors as used for the symbols). 



breakdown by instrument. We also present the cumula- 
tive histograms for various models with Cq included. In 
good seeing the bulk of the DIPSF distribution tends to- 
ward that of the DIMM, consistent with the Kolmogorov 
model. At the poor seeing end we see a migration of the 
DIPSF distribution toward lower FWHM values than the 
DIMM seeing. We find that a von Karman model with an 
outer scale value of Cq — 25 m bounds the distributions. 

4.1. Outliers and observer bias 

Some nights result in DIPSF FWHM values that are 
significantly higher than the corresponding DIMM and 
guide camera FWHM values. In the vast majority of 
cases these are MagIC data points, and all of these are 
traceable to nights on which the observer is known to use 
out of focus images in order to increase the total detected 
source flux and thus signal-to-noise ratio. The remainder 
of the offending points arise from engineering runs. In 
all of these cases we have removed the entire night's data 
from the analysis. The numbers presented above and in 
table [1] have these omissions taken into account. 

It is worth considering whether an additional bias is 
introduced by the observers themselves. Typically as- 
tronomers may drop an imaging program in favor of a 
backup spectroscopic program in the event of poor see- 
ing, and this may artificially bias our results in favor 
of nights with good seeing. To test the robustness of 
our results we compare the average DIMM seeing for 
our subset of nights with that for all DIMM data taken 
over the same period. The median DIMM seeing for 
our data set is 0'.'625 and that for the whole of the 
same time period (July 2007 - December 2008) is 0'.'675. 
T'hus we conclude that there is indeed an observer bias, 
and the median seeing reported above is slightly op- 
timistic for the year in question. However, when we 
compare to other periods, we find that the DIMM see- 
ing across our subsample is quite typical of long peri- 
ods at Las Campanas. For example, from September 



2005 - September 2008 the m edian DIMM seeing was 
0'.'62 ()Thomas-Os ip et al.ll20 1(]D: For the entire vea r 2007 
the DIMM seeing was 0"65 (|Tho mas-Osi p et all|2008( ); 
From April 200 5 - December 2006 the median DIMM 
seeing was 0'.'65 (|Thomas-Osip||2007f) . Thus we conclude 
that our estimate of Magellan image quality is robust for 
average conditions average, and that 2008 was an unusu- 
ally poor year. It is worth noting that 2008 was a La Nina 
year, and one of the wettest in the southern Atacama for 
a decade. 

4.2. Factors affecting the deviation of Magellan image 
quality from DIMM seeing 

It is clear from Fig. |4] that the Magellan imaging 
FWHM is consistently better than the DIMM seeing and 
the guide camera FWHM. In particular we find that the 
Magellan telescopes do much better than the DIMM see- 
ing when the seeing is poor. In order to provide a more 
powerful analysis, let us define a "seeing residual" that is 
the difference between the DIMM seeing and the science 
image quality, r^ = FWHMdimm - FWHMgcicnco, such 
that a positive r-g implies a better science image quality 
than the DIMM seeing. This gives us a direct measure of 
the improvement Magellan affords over the seeing for any 
given observations. We see no large-scale trend of see- 
ing residual with telescope pointing and conclude that 
there are no major directional impediments to the image 
quality obtained at Magellan. 

With the image quality database in place we are in 
a powerful position to begin exploring the effects of 
weather conditions on both the DIMM seeing and sci- 
ence image quality. It is known t hat the seeing corr elates 
weakly with strong wind speeds (|Consortiumll2006l ) . Fig- 
ure [H] illustrates the change of seeing residuals with wind 
speed. While there is always a weak tendency for r^ > 0, 
this becomes strongest at the highest wind speeds (and 
worst seeing conditions). In Figure [7] we illustrate the 
seeing residuals with wind vector in histogram form. No 
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Fig. 5. — Left: Cumulative distributions for the DIPSF, DIMM and guide camera FWHM measurements, von Karman model predictions 
are shown for two cases: £o = 25 m (dot-dashed) and Cq = 1000 m (triple-dot-dashed). Right: The instrumental breakdown of the DIPSF 
FWHM cumulative distribution (IMACS solid; LDSS3 dashed; MagIC dotted). To facilitate comparison with the left hand figure we show 
the DIPSF distribution for the complete dataset in grey. 



trend is seen with wind direction, beyond the fact that 
the prevaiUng wind direction at Las Campanas is North 
Easterly. Thus Magellan telescopes do relatively bet- 
ter than the DIMMs when the DIMM seeing is poor (as 
noted previously), and these situations often correspond 
to high wind speeds. But is this due to a genuine im- 
provement in the Magellan seeing, or just a worsening 
of the DIMM seeing? To help answer this we plot the 
seeing residual against both the DIPSF (blue triangles) 
and DIMM (red crosses) FWHM values in Fig. [5] We 
clearly see that the residuals are more strongly correlated 
with the DIMM values (Spearman's rank p = 0.69 with a 
probability p < 0.005 of occurring at random) than they 
are with the DIPSF values (p = 0.05; p = 0.18). Thus 
the relative improvement in Magellan image quality in 
poor seeing conditions is found to be, to first order, a 
local effect at the DIMM towers. The DIMM method 
is insensitive to tracking errors and wind shake, and we 
suggest that the high winds may cause local turbulence 
through interaction with the clamshell dome. The dome 
opens almost due south (bearing 177), and so will in- 
teract with the prevailing wind upstream of the DIMM 
apertures. 

We now remove all datapoints for which the wind speed 
exceeds 10 m s~^. This yields the cumulative image qual- 
ity distributions shown in Fig. [HI Removing the high 
wind speed data primarily improves the DIMM distri- 
bution, with a smaller improvement seen in the DIPSF 
distribution. The median DIPSF image quality drops 
from 0'.'575 to 0'.'553, whereas the median DIMM image 
quality drops from 0'.'625 to 0'.'594. The von Karman 
model distributions are defined by the DIMM distribu- 
tion and are therefore also affected. We see from the 
figure that the DIPSF distribution still cuts across a 
range of von Karman models. We note that an outer 
scale of £o ~ 25 m (consistent with other sites) gives 
a conservative estimate of the image FWHM at seeing 



values above 0'.'6. At low seeing values it seems likely 
that we are beginning to hit limitations imposed by the 
telescope optics. As discussed above (see foot of Sec- 
tion [3]), the LDSS3 and MagIC ports are shghtly better 
than the IMACS port, but in general we would expect 
even the best seeing conditions (~ 0'.'25) to return image 
FWHM of ~ O'.'a - 0'.'4 on the instrument ports, and this 
is observed in Fig. [9l 
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Fig. 6.— Seeing residual, r^ = FWHMdimm - FWHMgcience 
with wind speed. 



4.3. Negative DIMM bias 

As mentioned in section [2?3l we found our 2009 DIMM 
data to be affected by optical aberrations introduced 
by mis-coUimation and the resulting loss of ability to 
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= FWHMdimm - 
FWHMscioncc with wind direction in histogram format. The poor- 
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Fig. 8.— Seeing residual, Vs = FWHMdimm - FWHMgcionco 
against DIPSF FWHM (black/blue triangles) and DIMM seeing 
FWHM (grey/red crosses). The residuals show a strong correla- 
tion with the DIMM seeing (Spearman's rank p = 0.69 with a 
probability p < 0.005 of occurring at random) and a much weaker 
correlation with DIPSF FWHM (p = 0.05; p = 0.18). 



achieve a proper focus. Optical aberrations are ex- 
pecte d to cause both negat ive and positive bias up to a 
20% (jTokovinin et al.l [20071) . While positive biasing has 
been demonstrated experimentally b efore by the TMT 
site testing team (jWang et al.l |2007() , no reports have 
been published of a negative bias. 

The 2009 data set consists of 103 nights worth of 
data between January and October 2009. There are 
456 DIMM-DIPSF matches over 28 nights. The 2009 
dataset is therefore roughly half the size of the 2007-2008 
dataset (876 points over 66 nights). Figure fTOl compares 
the 2009 DIPSF and DIMM data seeing distributions 
to those used in the above analysis. We found that in 
2009 the DIMM profile was biased towards better seeing 
in good conditi ons and worse seeing in bad conditions. 
As discussed in iTokovinin et al.l ()2007D , a negative bias 
can occur for high altitude turbulence under certain con- 



ditions of optical aberration. In general at LCO, good 
seeing conditions occur when there is little ground layer 
seeing and the majority of the atmospheric turbulence 
occurs high up in the atmosphere. In this case, it is 
possible for the DIMMs to show a negative bias, depend- 
ing on the type of optical aberration present. Bad seeing 
conditions are more of a mix in terms of where the turbu- 
lence is located and thus may create more positive than 
negative biases. 

We used the MASS C^^{z) dz profiles and DIMM see- 
ing to explore the possibility of a correlation between the 
altitude of maximum turbulence and the seeing. Fig. [TT] 
plots the altitude of maximum turbulence as a function 
of seeing: the point size is an indication of how dominant 
the turbulence is in that layer (i.e. how much larger the 
maximum is than the secondary maximum in the profile) . 
We can see that for the dataset including the biased 2009 
data (right panel), the 16 km altitude bin clearly dom- 
inates in good seeing while this is less prominent in the 
unbiased dataset. Furthermore, there is also more tur- 
bulence seen in the 8 km bin in the biased dataset. Thus 
the negative bias detected is indeed associated with high 
altitude turbulence. 

4.4. Implications for observations in the NIR 

Co is independent of wavelength and thus the gain in 
image quality for a large telescope as we move redward 
is somewhat greater than we would expect for the strict 
case of Kolmogorov turbulence - see Fig. [121 If we adopt 
an average seeing of 0'.'65 and assuming £o = 25 m we 
would estimate that for a large telescope the ideal image 
quality (assuming perfect optics) at 500 nm is 0'.'52. In 
reality, we are degraded somewhat due to imperfect op- 
tics, as seen from the Shack-Hartmann system. If we ex- 
trapolate further, we estimate that the ideal image qual- 
ity for a large telescope should be 0'.'41 at 1/im, 0'.'29 at 
2^m, and 0'.'26 at 5/im. 

5. CONCLUSIONS 

We have presented a simple method for measuring the 
outer scale of turbulence that can be used at any obser- 
vatory site that has a large aperture telescope used regu- 
larly for imaging, and a DIMM device for measuring the 
seeing. The Las Campanas Observatory site offers truly 
excellent seeing conditions. The effective "seeing" for 
the Magellan telescopes is considerably underestimated 
by the DIMMs. While the DIMMs give a median seeing 
of 0'.'625, the median Magellan science image FWHM is 
0'.'575. We find the dominant cause of this difference 
is local to the DIMM: high wind velocities introduce 
local turbulence, perhaps through interaction with the 
clamshell dome. There is weak evidence for a directional 
infiuence on the DIMM image degradation at wind di- 
rections perpendicular to the orientation of the opening. 
However, removing points with wind speeds > 10 ms"-'^ 
we find that the Magellan image quality is still signifi- 
cantly better than the DIMM seeing. From the cumula- 
tive distributions of the DIPSF and DIMM FWHMs we 
estimate that the site has an average turbulence outer 
scale of Co ^ 25 m, similar to that measured at other 
sites. 

We were unable to include our 2009 data due to a mis- 
coUimation of the Cerro Manqui DIMM. This leads to 
both positive and negative biases in the DIMM data for 
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Fig. 9. — The cumulative image quality distributions for the DIPSF (solid black) and DIMM (dotted black) data for those data points 
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Fig. 10.— Cumulative FWHM distribution for 2009 for DIPSF 
(solid black line) and DIMM (dotted black line). The distribution 
for the data from July 2007 — December 2008 (as used in the bulk 
of this paper and plotted in Fig. [Sj is shown in grey for reference. 



the year, and is the first time that such negative biases 
have been recorded, although they are predicted. These 



negative biases occur when the upper layers of the atmo- 
sphere dominate the atmospheric turbulence profile. 
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